Universal scaling in transport out of equilibrium through a single quantum dot using 

the noncrossing approximation 
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The universal scaling behavior is studied for nonequilibrium transport through a quantum dot. 
To describe the dot we use the standard Anderson impurity model and use the non-equilibrium non- 
crossing approximation in the limit of infinite Coulomb repulsion. After solving de Hamiltonian, we 
calculate the conductance through the system as a function of temperature T and bias voltage V 
in the Kondo and in the mixed valent regime. We obtain a good scaling function in both regimes. 
In particular, in the mixed valent regime, we find excellent agreement with recent experiments and 
previous theoretical works. 



I. INTRODUCTION 

Scaling laws and critical phenomena are intertwined 
concepts through renormalization-group (RG) arguments 
as scaling ideas describe how physical quantities change 
under RG transformations. Scaling of physical quantities 
is particularly interesting because they reveal universal 
behaviour within certain global constraints. The power 
of universal scaling functions is that they associate seem- 
ingly different physical phenomena which shown similar 
behaviour when properly scaled with the aid of some pa- 
rameters. Specifically, zero-bias conductance through a 
single channel quantum dot (QD) coupled to a reservoir 
of conduction electrons as a function of temperature can 
be accurately fitted with an empirical formula extracted 
from the numerical renormalization-group calculations, 
Ge(T) |l|. This is due to the fact that at low tempera- 
tures and in the Kondo regime there is only one relevant 
energy scale, the Kondo temperature, Tr. Experimen- 
tal realizations of QD devices have reported a success- 
ful match between the obtained zero-bias conductance, 
appropriately scaled by TV, and the empirical formula 

g e (t) pig. 

When a finite bias voltage V is applied between the 
leads coupled to the QD, the system is no longer in equi- 
librium and a new energy scale emerges, eV (electron 
charge times voltage). It is, then, important to restore 
the corresponding universal scaling behavior. From the- 
oretical considerations on the single QD, it is expected 
that conductance as a function of temperature and bias 
voltage, G(T, V), may be described by a scaling relation. 
This scaling function has a quadratic temperature and 
voltage dependence, in the low range of T and eV 
compared to Tk- 

In a series of experiments in a spin- 1/2 Kondo dot, 
GaAs, Grobis et al. [(| measured the low temperature 
and low bias transport properties obtaining a quadratic 
power law in agreement with theoretical predictions. In 
addition to this, they suggested a universal function to 
describe the non-equilibrium conductance with the aid 



of only two scaling parameters even in the range of large 
temperatures as compared to the Kondo one. The scaling 
parameters are Tk and the height of the resonant peak 
in the limit that T — > 0, G . 

Recent experiments made by Scott et al. [3] on non- 
equilibrium transport through single molecule transistor 
in the Kondo regime have tested the universality of the 
scaling function of Grobis. The authors also report a 
quadratic form of the scaling exponents in the low energy 
range 

On the other hand, by using non-equilibrium renormal- 
ized perturbation theory (RPT), Rincon et al. [8| stud- 
ied the proposed universal law within the Kondo and the 
mixed valence regimes. The authors found appropriate 
scaling relation in both regimes when using the suggested 
scaling law. In addition, in the case of the mixed valence 
regime the parameters involved in the scaling formula 
match with the experimentally obtained for GaAs. 

In this work, we study the scaling properties of G(T, V) 
through a spin- 1/2 QD by using the non-equilibrium non- 
crossing approximation. We obtain the conductance as a 
function of temperature and bias voltage within the dif- 
ferent regimes of the model describing the QD. In what 
follows we discuss the observed scaling relation for each 
parameter range. The results for conductance and the 
extracted scaling parameters are compared with experi- 
ments and previous theoretical works [6rl8(. 



II. MODEL AND METHOD 

To describe the system in which the QD is coupled to 
two leads we use the spin 1/2— Anderson Hamiltonian: 



^2 E d n da + U n dt n di - 
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Here c\ va creates a conduction electron with momen- 
tum k and spin a in the lead v where v = L, R labels 
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the left and right leads. d\ creates an electron in the dot 
and rida is the number operator for a given spin in the 
QD. The energy of a single electron in the dot and the 
Coulomb repulsion are represented by Ed and U respec- 
tively. The hybridizations between the leads and the QD 
are given by Vk v and the lead-dot coupling strengths are 



A„( W )=7ry;o(w 



r), 



(2) 



which we consider to be step functions A u (u) = 
A l/ 0(D — \uj\) with an energy cut-off D of the order of 
several times A„. 

The bias voltage V is directly related with the differ- 
ence between the chemical potentials /i^ and [Ib. of the 
leads, eV = (1l — Mfl an d we set Ml = eVA^/A where 
A = A^ + Ar represents the total coupling between the 
leads and the QD. In this way, in the case of symmetric 
couplings, Al — Afl, the chemical potentials on both 
sides of the dot are related by ml = —p,R = eV/2. 

In a generic experimental configuration there are no 
restrictions on the relative couplings to the leads and 
the current through the dot can be expressed in terms 
of the lesser and greater Green functions of the QD Q . 
Otherwise, in the case of a proportional coupling, A^ = 
XAr, the current through the dot may be obtained solely 
in terms of the QD spectral function [Ic} 



I(V) = -AA J dujp d (oj) (/(w - ml) - /(w - Mfl)) , (3) 

where A = 4A^Ai{/A 2 represents the asymmetry fac- 
tor of the couplings between both leads and the dot. 
The function f(u) = 1/(1 + e^) is the Fermi dis- 
tribution and the spectral function Pd{^>) is given by 
Pd{u) = —ImG' dcr (uj)/ir. where G r da (oj) is the retarded 
Green function of the dot. 

In this work, we calculate pd(w) using the non crossing 
approximation (NCA) method in its strong coupling (SC) 
limit U — > oo. The NCA is a self-consistent conserving 
scheme in which the local operator d a is expressed by a 
combination of auxiliary boson and pseudo-Fermion par- 
ticles (Til. H3|. It yields an accurate quantitative descrip- 
tion of the equilibrium single channel Anderson model 
down to low temperatures, although it does not capture 
the correct Fermi-liquid behavior. As an additional ad- 
vantage of the method, the NCA can be gen eralized out 
of equilibrium in a straightforward way [Hi and it 
results in a reliable approximation for quantities involv- 
ing the Pd(w) such as the conductance as a function of 
temperature and bias voltage flil. [l5|. In the numerical 
procedure to solve the NCA equations, we have used a 
set of self-adjusting meshes (rather than the fixed one 
used in Ref. [l3j]) to describe the spectral densities and 
lesser Green's functions of the auxiliary particles. The 
procedure guarantees the resolution of the sets of inte- 
gral equations for the auxiliary particles self-energies to 
a high degree of accuracy. 



As we mentioned in the previous section, the purpose 
of this work is to test the experimentally observed scaling 
relation^ 0] and to provide a theoretical complement to 
the recent RPT study Q. 

The proposed universal scaling relation for the con- 
ductance as a function of temperature and bias voltage 
is given by the following expression 



G(T,V) 
G E (T) 



1 - 



acT (eV/ kTk) 2 



1 + ( 7 /q - l)c T (T/T K y 



(4) 



where Ge(T) is the empirical formula derived from a 
fit to numerical renormalization group conductance cal- 
culations 



Ge{T) = 



Go 



1 + (2 1 /* - 1) {T/T K f 



(5) 



and 



impurity. The constant 



s = 0.21 for a spin-^ 
ct = 5.49 is given by the low temperature expansion 
of Eqs. (0| and ([5]). The constants a and 7 represent the 
curvature at T — > and the evolution with increasing 
temperature of the resonant peak of the G(T, V), respec- 
tively. In a more general case when 7 = j(T) ~ l/T, 
the low temperature limit of the Ansatz given by Eq. (0]) 
agrees with the scaling Ansatz discussed by Hettler et at. 
[la | ( see Fig. 3b). This means that the proposal Eq. (0| 
in which 7 is not temperature-dependent represents the 
simplest approximation for "f(T). 

The empirical law, Ge(T), allows a definition of Tk as 
the temperature at which the equilibrium conductance is 
equal to half its zero temperature value, G(Tk) = Go/2. 

The scaling coefficients a and 7 are obtained and dis- 
cussed in the next section within different regimes of the 
Anderson Hamiltonian given in the Eq. ([T]). 



III. NUMERICAL RESULTS 

In the strong coupling limit, U —¥ 00 forbids double 
occupancy of the dot, therefore the number of electrons 
in the QD is restricted to be (n) < 1. In this limit there 
are tree different regimes of the QD, the Kondo- Regime 
(KR), the Mixed- Valence regime (MV) and the Empty- 
Orbital one (EO). In terms of the occupancy of the dot, 
the KR is characterized by (ukr) ~ 1- On the other 
hand, in the MV some charge fluctuations are allowed 
and therefore (umv) < 1- 

In the following we present the numerical results for the 
differential conductance as a function of temperature and 
bias voltage within the Kondo and mixed valence regimes 
in the U — > 00 limit. We focus on the case of symmetric 
coupling, A = 1, between the lead and the QD. This 
assumption is justified by the experimental setup used in 
the measurement of conductances for GaAS, that keep 
the two couplings nearly equal, [6|. 
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From now on, the total coupling A is taken as the unity 
of energy and we set the bandwidth to be D = 10A. 

The reader might wonder if a finite U expected in re- 
alistic systems affects our conclusions. We show that in 
the Kondo regime, to take infinite U is not an essential 
approximation. Below we will show that in the mixed 
valence regime, our results coincide with those of alter- 
native treatments developed for small U. In order to 
clarify the validity of the U infinite approximation con- 
cerning the scaling properties we compare in Fig.([T]) the 
equilibrium conductance G(T) /Go as a function of T/Tk, 
keeping constant the valence of the QD, for several values 
of the Coulomb repulsion U. To make a valid comparison, 
we use the finite (Ty| and infinite versions of the NCA well 
inside the Kondo regime with an occupation of the QD 
within the range 0.95 < (n) < 0.96. 

In the numerical evaluations we normalized the con- 
ductance to Go = G(Tq,0), where To is the lowest tem- 
perature that can be reached within our NCA calcula- 
tions (T - 0.005T K - 0.01T K ). 
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Figure 1: Equilibrium conductance versus T/Tk for several 
values of U. The valence of the QD is set around of 0.95. Tk 
is denned as the temperature such that G(Tk) = Go/2. 

From Fig.(TT]) it is clear that within the Kondo regime, 
in which the charge fluctuations are frozen, the only rel- 
evant scale is given by the Kondo temperature, Tk- In 
Ref. [l?} , Gerace et al. performs a detailed calculations of 
scaling behaviour with finite U values within NCA. 

As we mentioned in the previous section, in what fol- 
lows, we obtain the conductance through the QD as 
a function of temperature and bias voltage within the 
strong coupling limit U — > oo in the Kondo regime and 
mixed valence one. 

Within the NCA these regimes are schematically indi- 
cated in the inset of the FigJ5] As it is shown the cases 
Ed < — 3 A and — 2 A < Ed < can be considered as 
describing the KR and MV regimes, respectively. 

In order to obtain the differential conductance G(V) = 
dl (V) I dV , we perform a numerical differentiation of the 
current in Eq. ((3]) for a given T. Having calculated Tk 
from the empirical law Ge (T) , the value of a in the uni- 



versal scaling function, Eq. ([J}, can be determined from 
its low temperature expansion 



G(T , V)/G ~ 1 - ac T (eV/nT K y 



(6) 



Notice that from RPT calculations in the strong cou- 
pling limit U — > oo, the value of the a parameter does 
not depend on the asymmetry A. This follows from Eq. 
(|7|), in which the strong coupling limit is given by the 
renormalized parameter u — U / (nA) — > 1 (l8| . 



g_ _ 7T 2 (i + 2u 2 ) /fcry 

Go 3 \a) 

4 - 3A + (2 + 3A)u 2 _ /eV 

X 



(7) 



Once a is obtained, the 7 parameter is found by plot- 
ting the scaled conductance (1 - G(T, V) /G(T, 0))/a„ 
versus (eV/ kTk) 2 where 



acT 



1 + (7/a - l)c T {T/T K y 



(8) 



A. Kondo regime 

We begin analyzing the case in which Ed = — 3A and 
the QD is symmetric connected to the leads, A^ = Ar. 
We obtain Tjc by a fit the NCA linear response conduc- 
tance (V — 0) with the empirical law Ge{T). In FigJ^l 
we compare the empirical law with the NCA result. The 
obtained conductance follows the empirical law very well 
in the whole range of temperatures and fixes the value of 
the Kondo temperature as Tk = 0.038A. 

When a finite bias V is applied the differential conduc- 
tance as a function of the bias voltage exhibit the usual 
Kondo peak centered at zero bias and it broadens fast 
as the temperature increases. The obtained G(T, V) are 
shown in Figj3]as a function of eV/ kTk for several tem- 
peratures below the Kondo one. It is noticed that in the 
limit of T — >• the conductance saturates at the usual 
value of 2e 2 /h for the Kondo regime. 

From Eq. ^ we extracted the value of a at very low 
temperatures as compared to Tk ( O.OITr- < T < 0.05Tr- 
) and over the range eV/ kTk < 0.4 as in the GaAs exper- 
iment. We obtain an average value a = 0.16 ±0.01. This 



value of a is consistent with the 



,,SAM _ 
Y SC 



3/(2tt 2 ) = 0.152 



predicted for the Anderson model in its strong U — > 00 
limit The inset of Figl3] shown the fitted conductance 
for T = 0.05T K - 

The values of 7 for each temperature are obtained from 
Eq[8] and are listed in Table fl] From the temperature 
dispersion, Table HI the average value of 7 is 7 ~ 0.9 ± 
0.2. The obtained values for a and 7 are major than the 
corresponding ones to the experiment. 
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T/T K 

Figure 2: Equilibrium conductance versus empirical law 
Ge(T) as a function of temperature within the KR regime, 
Ed = — 3A. The squares represent the NCA conductance 
while the solid line is for the empirical law defining Tk = 
0.038A. In the inset the different regimes are schematically 
indicated as a function of Ed/ A- 
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Figure 3: Differential conductance as a function of eV/nTpc 
for several temperatures below Tk and for Ed — — 3A. The 
inset shows the linear regression curve that fitted G(T, V) very 
well up to (eV/KT K ) 2 ~ 0.2 at T = 0.05T K . The obtained 
value for a is a = 0.16 ± 0.01 



Table I: Shown are the calculated values for a v and 7 as a 
function of temperature for Ed = — 3A. Temperature in units 
of T K . 



Temperature a v 



7 



0.05 
0.10 
0.20 
0.30 
0.40 
0.50 



0.743 1.33 
0.643 0.93 
0.402 0.87 
0.270 0.77 
0.178 0.77 
0.132 0.73 



viation of the Ansatz scaling due to the temperature dis- 
persion of 7 was also found by using RPT (Ref. [§[). Be- 
sides, our results agree with the obtained by using RPT, 
a = 0.151 and 7 = 0.8, given in Ref. [|[. 
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Figure 4: Scaled conductance (1 - G(T, V) /G(T, 0))/a v as 
a function of (eV/ kTk) 2 for several temperatures in the case 
of Ed = — 3A. The temperatures are given in units of Tk- A 
collapse of all curves into a linear y = x is obtained. 



B. Mixed valence regime 



In FigQJ we plot the scaled conductance (1 — 
G(T, V)/G(T, 0))/a v as a function of {eV/nT K ) 2 for sev- 
eral temperatures using the average calculated values of 
a and 7. The collapse of all curves in a straight line with 
slope 1 at low voltages might suggest that the Ansatz 
scaling is appropriate. However, an important point to 
note is that this collapse require different values of 7 for 
different temperatures. The dependence of 7 with the 
temperature could indicates the breakdown of the scal- 
ing Ansatz given by Eq. Q in the Kondo regime. The 
decrease of 7 with increasing temperature, as was men- 
tioned in the previous section, agrees with a more general 
relation in which 7 is proportional to 1/T (l3| . This de- 



We now focus on the nonequilibrium transport through 
the dot for E4 = —2 A. In this regime of the parameters 
the occupancy in the QD is (n) ~ 0.8 (inset of Fig(2]) 
so some charge fluctuations are allowed. From the con- 
ductance at zero bias and its fitting with the Ge(T), 
the corresponding Kondo temperature is found to be 
T K = 0.158A. As it is shown in FigE] the fit of the 
NCA conductance with the empirical law is very good in 
the whole range of temperature. 

Figure [S] shows the temperature evolution of the con- 
ductance as a function of the applied bias. As it can 
be seen from a comparison between FigsJ3] and [HI the 
curvature of the conductance at T — > is smoother in 
the present case than for Ed = —3 A. It is noticed that 
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Figure 5: Equilibrium conductance versus empirical law 
Ge(T) as a function of temperature within the MV regime, 
Ed = — 2A. The squares represent the NCA conductance 
while the solid line is for the empirical law defining Tk = 
0.158A. 



the conductance saturates at a value less than the corre- 
sponding to the Kondo regime (2e 2 /h) due to the allowed 
charge fluctuations. 

Following the procedure already described in the pre- 
vious section, from the G(T, V) at T = T we obtain 
a = 0.11±0.01. This value of a is in a surprising agree- 
ment with the one obtained by Grobis et al. in the GaAs 
experiment [6|. Moreover, this is qualitatively the same 
feature that was found within the RPT [§|. 
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Figure 6: Differential conductance as a function of eV/ kTk 
for several temperatures below Tk for Ed = — 2A. The inset 
shows the linear regression curve that fitted G(T, V) very well 
up to (eV/nT K ) 2 ~ 0.3 at T = 0.05T K . The obtained value 
for a is a = 0.11 ±0.01. 



and finite U in contrast to our case. 



Table II: Shown are the calculated values for a v and 7 as a 
function of temperature for Ed = — 2A. Temperature in units 
of T K . 



Temperature a v 



7 



0.05 
0.10 
0.20 
0.30 
0.40 
0.50 



0.575 0.514 

0.505 0.502 

0.340 0.498 

0.225 0.485 

0.158 0.462 

0.115 0.452 



Fig[7] shows the scaled conductance (1 — 
G(T,V)/G(T,0))/a v versus {eV/nT K ) 2 for several 
temperatures using the average calculated values of 
a and 7. We notice that the Ansatz for the uni- 
versal scaling relation works very well up to at least 
{eV/nT K ) 2 < 0.5 and for 0.05 < T/T K < 0.5. As in the 
case of KR it can be seen that all curves collapse into 
the linear y = x. 
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Figure 7: Scaled conductance (1 - G(T, V)/G(T,0))/a v as a 
function of (eV '/ kTk) 2 for several temperatures in the case of 
Ed ~ — 2A. The temperatures are given in units of Tk- A 
collapse of all curves into a linear y = x is obtained. 



IV. CONCLUSIONS 



The 7 values within this regime, and more important 
their dispersion are, also reduced in comparison with the 
respectives obtained for the KR. Table [TT] shows the val- 
ues of a v and 7 from the fit of conductance for each tem- 
perature. The average value of 7 from the temperature 
dispersion is 7 ~ 0.48 ±0.03, in excellent agreement with 
the experimental value 7 ~ 0.5±0.1, @. This value agree 
with that found by Rincon et al. in this regime although 
they used a symmetric Anderson model, Ed = —U/2, 



In the previous sections we analyzed the Ansatz of Gro- 
bis et al. for the universal scaling of the conductance 
out of equilibrium through a quantum dot described for 
a spin 1/2— Anderson model. We performed a non- 
equilibrium non-crossing approximation study of the con- 
ductance within the Kondo and Mixed Valence regimes. 
In both regimes we obtain an appropriate scaling rela- 
tions and our results agree with the previous ones by us- 
ing the non-equilibrium renormalized perturbation the- 
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ory. However, in the Kondo regime, different values of 7 
for different temperatures T should be used. This sug- 
gests that the fitting formula is inappropriate, at least in 
the Kondo regime, and the dependence of 7 with tem- 
perature should be considered. 

This is the same observation obtained by Rincon et al. 
[1| even the method and the model are quite different that 
the ones used in this work. Both techniques, non-crossing 
approximation and renormalized perturbation theory, are 
suitable methods for study the transport properties out 
of equilibrium for the Anderson model. The agreement 
represents a good benchmark between both of them. 

The quasi-constant value of 7 in the mixed valence 
regime indicates that the Ansatz given by Eq. |4]) is an 
appropriate scaling relation in this regime. The scaling 
parameters depend on the regime of the dot. In partic- 
ular, if valence fluctuations are allowed, the parameters 



are in excellent match with the corresponding ones ob- 
tained in experiments on GaAs QD's. Furthermore, this 
agreement with the experimental values of the a and 7 
parameters suggests that the system in the experiments 
of Grobis et al. presents some valence fluctuations char- 
acteristic of the mixed valence regime. 
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